Pseudotensors and quasilocal energy-momentum 
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Early energy-momentum investigations for gravitating systems gave reference frame dependent 
pseudotensors; later the quasilocal idea was developed. Quasilocal energy-momentum can be deter- 
mined by the Hamiltonian boundary term, which also identifies the variables to be held fixed on the 
boundary. We show that a pseudotensor corresponds to a Hamiltonian boundary term. Hence they 
are quasilocal and acceptable; each is the energy-momentum density for a definite physical situation 
with certain boundary conditions. These conditions are identified for well-known pseudotensors. 



On ' 
On : 
ON ■ 
t— i ■ 

5h ! 

<: 
<n : 

(N : 
> ■ 
o : 
. 

o ■ 

On ■ 

o : 

00 . 
On ■ 

"o ! 
i ' 
bo 



Energy-momentum can be regarded as the most funda- 
mental conserved quantity, being associated with a sym- 
metry of the space-time geometry. From Noether's theo- 
rem and translation invariance one obtains a tensor, T MI/ , 
which describes the density of energy- momentum and de- 
fines a conserved energy-momentum as a consequence of 
its satisfying the differential conservation law d u T^ v = 0. 
But it is not unique, since one can add an arbitrary 
"curl" — which shifts the zero of energy- momentum. The 
gravitational response, however, via Einstein's equation 
G^ v — kT mi/ , detects the total energy- momentum den- 
sity for matter and interaction fields, and thereby re- 
moves the classical uncertainty in the energy-momentum 
expressions — the gravitational field is trivial only if the 
total energy momentum vanishes. Thus it seems ironic 
that the proper identification of the contribution to the 
total energy-momentum from the gravitational field it- 
self has proved to be so elusive. It was natural to ex- 
pect that, since these other sources exchange energy- 
momentum with the gravitational field locally, gravity 
should also have its own local energy-momentum density. 

Attempts at identifying an energy-momentum density 
for gravity, however, led only to various energy momen- 
tum complexes which are pseudotensors, including those 
of Einstein Gj, Papapetrou Bergmann ||, Landau 
and Lifshitz || , M0ller || , and Weinberg Q . Pseudoten- 
sors are not covariant objects; they inherently depend on 
the reference frame, and thus by their very nature can- 
not provide a truly physical local gravitational energy- 
momentum density. Indeed any such quantity is pre- 
cluded by the equivalence principle itself, since a gravita- 
tional field should not be detectable at a point. Conse- 
quently many have criticised the whole idea; e.g., an in- 
fluential textbook states: Anyone who looks for a magic 



formula for "local gravitational energy-momentum" is 
looking for the right answer to the wrong question. Un- 
happily, enormous time and effort were devoted in the 
past to trying to "answer this question" before investi- 
gators realized the futility of the enterprise. |7j p. 467. 
Hence the pseudotensor approach has been largely aban- 
doned (although interest continues, see e.g., Q). A new 
idea, quasilocal (i.e., associated with a closed 2-surface), 
was proposed and has become widely accepted. In view 
of the role of gravity in identifying energy, this then be- 
comes the most fundamental notion of classical energy- 
momentum. There have recently been many quasilocal 
proposals [9 10 1. Various criteria have been advanced 
(see, e.g., [ll|), in particular good limits to flat space- 
time, weak field, spatial infinity (ADM), and null infinity 
(Bondi). However it has now been recognized that there 
are an infinite number of expressions satisfying these re- 
quirements jn| . Clearly, additional principles and crite- 
ria are needed. The Hamiltonian formalism includes such 
a principle, with the Hamiltonian boundary term deter- 
mining both the quasilocal quantities and the boundary 
conditions |l^,^3|. Here we show that this Hamiltonian 
approach to quasilocal energy-momentum rehabilitates 
the pseudotensors. 

We identify energy with the value of the Hamiltonian. 
The Hamiltonian for a finite region, 



H(N) 



N"n L 



B(N), 



(1) 



S=<9£ 



generates the spacetime displacement of a finite spacelike 
hypersurface S along a vector field N^; it includes a sur- 
face and a boundary term. Noether's theorems guarantee 
that Tifj, is proportional to the field equations. Conse- 
quently the value depends only on the boundary term 
B, which gives the quasilocal energy-momentum. But the 
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boundary term can be modified. (This is a particular case 
of the usual Noether conserved current non-uniqueness.) 
Indeed it is necessary to adjust B to give the correct 
asymptotic values Fortunately, B is not arbitrary. 

A further principle of the formalism controls its form: 
choose the Hamiltonian boundary term B so that the 
boundary term in SH vanishes when the desired fields are 
held fixed on the boundary (technically necessary for the 
variational derivatives to be well defined). Hence, we find 
a nice division: the Hamiltonian density Ti.^ determines 
the evolution and the constraint equations, the bound- 
ary term B determines the boundary conditions and the 
quasilocal energy-momentum. There still remain many 
possible boundary condition choices Consequently 
there are various kinds of energy, each corresponding to 
a different choice of boundary condition; this situation 
can be compared with thermodynamics with its various 
energies: internal, enthalpy, Gibbs, and Helmholtz. 

For geometric gravity theories including general rel- 
ativity, the Hamiltonian can be succinctly obtained in 
terms of differential forms |l7|| , which readily displays 
the boundary term and its connection to pseudoten- 
sors. Since we wish to connect with traditional work 
we shall use a coordinate (holonomic) basis dx a and 
sometimes the (Hodge) dual basis n a " := *(dx a A • • •). 
The simplest analysis uses the connection along with 
the metric as the dynamic variables. The curvature 2- 
form is given by £l a p := du> a p + cl> q 7 A ui J p, where 
uo a p — r Q p-ydx 7 is the connection one-form. The La- 
grangian density is the Einstein-Hilbert scalar curvature 
4-form C = fl a p A nj 3 = R^/=g~d A x. The Hamilto- 
nian 3-form can be constructed according to the pat- 
tern L = qp — H by contracting the Lagrangian 4-form 
with the time evolution vector field (dropping indices for 
brevity): 

i N C = i N n A 77 + fi A i N r/ = £ n uj A n - H(N); (2) 

here the time derivative is given by the Lie derivative (on 
form components £n :— dipj + ijvrf); the Hamiltonian 3- 
form (without discarding total derivatives) is 



H(N) = — fi A i N r, - i N uDr) + dB M (N). 



(3) 



The Hamiltonian density term includes N^TC^ = —£l a pA 
N^-qa 13 ^ = 2N^G V ^r\ v , a covariant expression which 
projects to the usual ADM Hamiltonian density (see, 
e.g., Ch 21 in 0), along with a frame gauge trans- 
formation generating term, i^coDr], which vanishes be- 
cause the connection is symmetric and metric compat- 
ible. When integrated over a finite spatial hypersur- 
face E, the value of the Hamiltonian comes from the to- 
tal differential term which (via the generalized Stokes 
theorem) yields a boundary term with the 2-form in- 
tegrand B M (N) = i N u a pr)J = N X M^» (l/2)dS M „, 
where dS^ := (l/2)e fiva pdx a A dx^\ the coefficients 



turn out to be the superpotential whose divergence gives 
the M0ller pseudotensor M . (Note the extreme directness 
of this derivation of M0ller's expression.) The variation 
of the Hamiltonian density (0) yields (with N fixed) 



SHm = field equation terms + dijv(&j Q p A -qj' 



(5) 



showing that the boundary condition implicit in (||) is of 
the Neumann type: the connection r ~ dg is to be held 
fixed. This calculation also reveals a serious deficiency: 
the boundary term in the variation of the M0llcr Hamil- 
tonian will not vanish with the standard asymptotics: 
6g~0(l/r), 8Y~0{l/r 2 ). 

An almost obvious alternative (in view of the problem 
just mentioned) is simply to replace the M0ller Hamil- 
tonian boundary term with Be '■= ^ a /3 A iNVa^ 3 — 
N x U x ^(l/2)dS^. Now the coefficients 



= (-grhxrd^ l-g(9» T g^ - gTsD], (6) 



are the Freud superpotential, whose divergence gives the 
Einstein pseudotensor [jl 18|. The Hamiltonian variation 
now contains a boundary term of the form di nC where 



C 



-u a fs A 5t,J = -T a Pl 8{^g-gP a )8Z{d*xU (7) 



= 2y 



'99 



^T a \ jx = 2^g^g^d w g^ ]x , (4) 



showing that this choice of Hamiltonian boundary term 
is associated with a Dirichlet type boundary condition: 
the contravariant metric density is to be held fixed. 

The possible forms of the Hamiltonian boundary term 
have been considered in some detail elsewhere Jl3| . There 
are various choices involved including: (1) the represen- 
tation or dynamic variables, such as the metric, orthonor- 
mal frame, connection, spinors; (2) the control mode: 
such as a Dirichlet or Neumann boundary condition for 
each dynamic variable; (3) the reference configuration: 
this is given by the field values that give vanishing energy- 
momentum (the standard choice is Minkowski space). (4) 
the displacement vector field N: which timelike displace- 
ment gives the energy? which spatial displacement gives 
the linear momentum? 

We have shown how two famous pseudotensors natu- 
rally arise from Hamiltonian boundary terms and how 
they are associated with field boundary values. In 
like manner we shall now explain how the Hamilto- 
nian boundary term approach to quasilocal energy- 
momentum rehabilitates all of the other pseudotensors 
so they can be recognized as legitimate. There is a di- 
rect relationship between pseudotensors and quasilocal 
expressions. Every pseudotensor corresponds to some ac- 
ceptable choice of boundary expression. Conversely, ev- 
ery boundary expression defines a pseudotensor Jl9| . 

We consider the pseudotensor idea in some detail: a 
suitable superpotential H fl uX = H^ vX ^ is selected and 
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used to split the Einstein tensor thereby defining the as- 
sociated gravitational energy-momentum pseudotensor: 

Ky/=gNH„ v := -N»y/=gGS + ~3 A (iV M ff/ A ), (8) 

where K — 8irGc~ and we have inserted a vector field to 
make the calculation more nearly covariant. The usual 
formulation is recovered by taking the components of the 
vector field to be constant in the present reference frame; 
then Einstein's equation, — kT^", can be rearranged 
into a form where the source is the total effective energy- 
momentum pseudotensor 

d x H^ x = 2 K (-#V := 2k(- S )*(V +T/). (9) 

An immediate consequence of the antisymmetry of 
the superpotential is that is a conserved current: 
d u [{—g) 1 / 2 T tl v ] = 0, which integrates to give a conserved 
energy-momentum, N^P^ := / N< J 'T^ v (-g) 1 ^(d s x) v . 
This should be contrasted with the covariant formula 

v„ v = W - r A , M r A " + t v Xv t v x = o, (10) 

which does not lead to a conserved energy-momentum 
unless r = (flat space). 

A minor variation on the preceding analysis results 
from choosing a superpotential with a contravariant in- 
dex: H^ vX = H^ vX \ A further variation: 

H" va := dpH^, (11) 

along with the symmetries H» av0 = H"^ a = 
and = o, leads to a symmetric pseudotensor — 

which then allows for a simple definition of angular mo- 
mentum, see @] §20.2. We can cover these options sim- 
ply by introducing a vector field, then, in the follow- 
ing computations, we can easily make modifications like 
N^H^ X — ► N^H^ uX . 

The energy-momentum within a finite region 

- P(N) :=- [ N<*T/^g-(d 3 x) v = 

/ [N"V=gfiG/ - V) - ±-d x (N»H^ x )} (d 3 x) v 

= N^+i B(N)=H(N), (12) 

JT, JS=dT, 

is seen to be just the value of the Hamiltonian. Note 
that Tifj, is the covariant form of the ADM Hamiltonian 
density, which has a vanishing numerical value, so that 
the value of the Hamiltonian is determined purely by 
the boundary term B{N) = -N^(l/2K)H^ x (l/2)dS v \. 
Thus for any pseudotensor the associated superpotential 
is naturally a Hamiltonian boundary term. Moreover 
the energy-momentum defined by such a pseudotensor 
does not really depend on the local value of the reference 
frame, it is actually quasilocal — it depends (through the 



superpotential) on the values of the reference frame (and 
the fields) only on the boundary of a region. 

Even more important, the Hamiltonian approach en- 
dows these values with a physical significance. To un- 
derstand the physical meaning of the quasilocalization, 
calculate the boundary term in the Hamiltonian varia- 
tion: 

- -L [ST a fl xN' t V=ggf , "S% r x tl + 6(N»H^)] dS Tp . (13) 

For the Einstein pseudotensor, we use the Freud super- 
potential (^) as the Hamiltonian boundary term in p2] ) . 
Then the boundary term in the Hamiltonian variation 
H has the integrand 5{y/^gg^ a 'N»)T a pxS^*, which 
shows not only that \/—gg® a is to be held fixed on the 
boundary, but also that the appropriate displacement 
vector field is N 1 - 1 = constant, and the reference con- 
figuration here (as well as in the other cases below) is 
Minkowski space with a Cartesian reference frame. 

This calculation is easily adapted to some other 
cases just by adjusting N* 1 . The Bergmann pseu- 
dotensor [|, given by 2 Ky ^T^ := d x (g^U 7 " x ), 
leads to the Hamiltonian boundary variation term 
SiV^g^g^N^pxd^dSrp, revealing that we are 
to fix V — 35^ 5 011 the boundary and use the dis- 
placement (co)vector = N® = constant. This last 
statement takes on a more proper geometric form in 
terms of an auxiliary (background) metric g^ having 
the Minkowski values in this coordinate system, then the 
desired displacement vector is N a = g aj N^. 

The Landau- Lifshitz pseudotensor H is slightly 
more complicated, being given by a weighted density 
2K{-g)T£Z := a A (V 3 55 P7 ^7 l/A )- The easiest way to 
handle this is to introduce, where necessary to obtain 
the proper geometric density weights, "extra" factors of 
(— g) 1 / 2 , the Jacobian factor for the flat metric (numer- 
ically constant in Cartesian coordinates). This leads to 
the conclusion that the displacement vector should be 
N» = guv (g/g) 1 / 2 N° and the quantity to be held fixed 
on the boundary is {—g)g^ rT g^^ ■ 

The three pseudotensors just discussed are associated 
with similar but distinct Dirichlet type boundary con- 
ditions which are algebraic in terms of the metric. On 
the other hand, the M0ller pseudotensor has a simple 
Neumann type condition. While the detailed physical 
significance of these conditions has not yet been probed, 
such an investigation seems straightforward. In contrast, 
the remaining pseudotensors in our survey are associated 
with more complicated boundary conditions. 

In the context of Eq. (pr|), Goldberg (2ll discussed the 
general form H^ a » p = H^ v H af} ~~H av including var- 
ious weighted densities. Because of the symmetries, the 
associated pseudotensor, <J^g~T^ v {H) :— d a dj3H^ avf3 , is 
guaranteed to be symmetric and conserved for all sym- 
metric H^ v . This pattern can nicely accommodate the 
Landau-Lifshitz version with H^ v := \J—gg^ v - More 
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generally, we note that r{H 1 +H 2 )-T{H 1 )- T(H 2 ) 
is also identically conserved leading to the more gen- 
eral pattern H» av P = H?" ' Hf - H[ lf} H% a + H£* ' Hf - 
Hf H^ a . With Hf — g a>3 we can now accommodate 
the pseudotensors of Papapetrou ||] (H^ u = yj—gg^") 
and Weinberg § (#f = -W + \g^ v h x x , where := 
9nu — 9fj,v and indices are raised with g^). The Hamil- 
tonian variations then lead to boundary conditions in- 
volving rather complicated combinations of 5T a ^ v and 
S{y/—gg^ v ) or Sh^, respectively. 

In summary, because of the very nature of the gravita- 
tional interaction and its elusive contribution, the local- 
ization of total energy-momentum has remained an out- 
standing fundamental puzzle. The earlier pseudotensor 
approach was considered to be unsatisfactory. A newer 
idea is quasilocal. Quasilocal energy-momentum can be 
obtained from the Hamiltonian. For a finite region it in- 
cludes a boundary term which plays the key role, deter- 
mining both the boundary conditions and the quasilocal 
values. Consequently there are (as in thermodynamics) 
many different physical kinds of energy, each correspond- 
ing to a different boundary condition. We have shown 
that every energy-momentum pseudotensor is associated 
with a legitimate Hamiltonian boundary term. Hence the 
pseudotensors are quasilocal and acceptable. Each is the 
energy-momentum density for some definite physical sit- 
uation. Via the Hamiltonian approach one can identify 
the necessary boundary conditions and thereby appre- 
ciate the physical significance of the associated energy- 
momentum quasilocalization. 

Our analysis reclaims the pseudotensor work of the 
past to its rightful place: concerning a special class of 
quasilocal energy-momentum. Moreover it is additional 
evidence that the Hamiltonian boundary term approach 
to energy provides an effective ordering principle for the 
various quasilocal energy- momentum expressions. 
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